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Abstract  This  data  report  contains  results  from  selected  time  series  and  spec¬ 
tral  analyses  of  the  turbulent  wind  and  temperature  measurements  performed  by 
Riso  National  Laboratory  during  the  AMADEUS  **5fmoke  and  Diffusion  Tests” 
of  Project  WIND,  Phase  IV.  which  took  place a  cooperative  research  oriented 
study  between  the  U.S.  Army  Atmospheric  Sciences  Laboratory  (ASL)  and  U  S. 
Department  of  Agricultural  Forest  Service  (USDAFS)  in  the  Meadow  Brook  Val¬ 
ley  near  Read  Bluff.  California,  during  the  period  between  3  September  and  7 
October  1987.  /..**£  - 

In  Part  I  of  this  study  (Riso-  M-2718.  January  1989).  Riso  National  Labora¬ 
tory  reported  sonic-anemometer  measurements  of  10-min  averaged  surface-layer 
scaling  parameters  such  as  surface  heat  flux,  shear  stress,  turbulence  levels  and 
atmospheric  stability  measured  at  two  locations  in  the  Meadow  Brook  Valley 
floor  accompanied  by  aerial  photography  of  the  valley-floor  smoke  puff  and  plume 
spread. 

The  present  study  (Riso-M-2861)  provides  time  series  plots  of  the  turbulent 
( 10  Hz  block-averaged)  wind  and  temperature  signals  as  recorded  by  sonic-anemo¬ 
meters /thermometers  at  the  7-m  level  above  the  Meadow  Brook  Valley  floor  during 
the  AMADEUS  trials.  The  time  series  are  further  processed  into  energy  spectra 
for  the  three  wind  components  (u',i/.u/)  and  fluctuating  temperature  (T')  and 
here  presented  together  with  their  relevant  scaling  parameters  calculated  by  the 
correlation  method. 

The  time  series  and  spectra  provide  flow  and  diffusion  modellers  of  the  AMADEUS 
experiments  with  an  insight  in  the  turbulent  scales  and  energies  most  responsi¬ 
ble  for  the  observed  flow  and  diffusion  processes.  Furthermore  they  provide  high- 
resolution  boundary-layer  flow  and  turbulence  measurements  for  model  simulation 
of  the  individual  experiments. 

All  data  have  been  transferred  to  ALS  on  IBM  PC-compatible  diskettes 


The  work  was  performed  under  Contract  No  DAJ A45-88-M-0325  with  the 
U  S.  Army  Research.  Development  and  Standardization  Group  -  U  K  under  su¬ 
pervision  of  Chief.  Environmental  Sciences  Branch,  Mr  J.C.  Comati. 
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1  Introduction 

1.1  Measurements 

& 

Besides  on-line  computer-system  calculating  the  10  min  mean  statistics  of  shear 
stress,  hej  flux,  variances  etc.  based  on  the  four  sonic  signals  as  described  in  part  l. 
additional  (high-frequency)  analog  recording  of  these  signals  were  taken  from  the 
central-valley  located  "sonic  B”  during  a  total  number  of  14  diffusion  tests,  with 
the  purpose  of  performing  time  series  and  spectral  analyses  of  the  surface-layer 
turbulence  quantities. 


Table  1.  .4  MADE  US  1981  sonic  spectra 


i  Speed:  2.38  cm/s  (15/16 
ch  1  ch  2  ch  3 

In)  (t>)  (u>) 

10  m/s  10  m/s  2  m/s 

(max  8  hrs) 
ch  5 
(T) 

5°C  /Volt 

AMADEUS  1987  Sonic 

Spectra 

FM  -  tape  j 

Date 

Start 

Stop 

Duration 

Spectra 

*  i 

23  Sep 

11:00 

15:57 

4  hrs  57  mm 

14:00  ■  15:32 

#  2 

24  Sep 

Wr- ^  ajj 

05:07 

8  hrs  08  min 

23:45  -  01:17  1 

#  4 

26  Sep 

13:32 

3  hrs  02  min 

12:00  •  13:32 

#  5 

27  Sep 

06:15 

12:15 

06:15  -  07:47 

#  6 

28  Sep 

14:52 

4  hrs  32  min 

10:20  -  11:52  ! 

#  ~ 

30  Sep 

06:22 

8  hrs  08  min 

06:32  -  08:04 

#  8 

30  Sep 

17:35 

01:43 

8  hrs  08  min 

19:25  -  20:57  ■ 

#  9 

01  Oct 

06:32 

14:42 

3  hrs  10  mm 

06:32  -  08:01 

#  10 

01  Oct 

06:32 

14:42 

8  hrs  10  min 

13:10  -  14:42 

#  11 

02  Oct 

06:40 

13:20 

06:40  -  08:12 

*  12 

02  Oct 

06:10 

8  hrs  40  min 

11:48  -  13:20 

#  13 

03  Oct 

06:40 

12:34 

5  hrs  54  min 

06.40  -  08-12 

4  14 

03  Oct 

14:10 

22:18 

19:30  -  2102 
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2  Data  Analysis 

2.1  Selection  of  runs 

An  overview  of  the  14  diffusion  tests  is  given  in  table  1. 


2.2  Digitizing,  block  averaging  and  time  series 
plots 

In  the  first  step,  the  four  raw  sonic  signals  from  each  experiment  were  digitized  at 
our  laboratory  at  a  rat  eof  10  Hz  and  calibrated  into  physical  units  of  wind  speeds 
U.  V',  and  temperature  T. 

Secondly,  each  raw  data  series  was  examined  and  quality  assured  before  a  se¬ 
quence  of  approximately  92  minutes  (55.000  scans)  was  selected  for  further  anal¬ 
ysis.  comprising  the  following  steps. 

1.  Calculation  of  the  V  :  V  and  W  mean  values  followed  by  a  rotation,  first 
around  the  y-axis.  next  around  the  c-axis  of  the  anemometer,  in  order  to 
obtain  externally  aligned  data  series  in  which  the  mean  values  (V)  and  (}V) 
are  equal  to  zero. 

2.  Calculation  of  the  symmetric  covariance  matrix  of  fluctuating  quantities: 


(tin)  (uv) 

(uw) 

(ut) 

M 

(»«■> 

(vt) 

(tan) 

M 

(«> 

based  on  the  entire  92-mi ri  aligned  subset  of  approximately  55.00U  scans 
where  n  =  V  -  (V),  v  =  V'  -  (V),  w  =  IV*  -  (W)  and  t  =  T  -  (T). 

3.  High-resolution  time  series  plots  were  made  of  the  aligned  three-wind  com¬ 
ponents  and  temperature  signal  for  inspection. 


2.3  Spectral  analysis:  velocity  and  temperature 
spectra 

From  each  experiment,  consisting  of  some  55.000  measurements  of  V,  V.  U'  and  T 
five  to  six  segments  of  equal  length  of  213  (8192  scans,  euqal  to  819.2  sec)  data 
points  were  selected.  These  were  in  turn  individually  Fourier  transformed  into 
•4094  spectral  power  values,  ranging  in  frequency  approximately  3.6  decades  from 
1/819.2  Hz,  2/819.2  Hz.  3/819.2  Hz  4/  ...  up  to  the  5  Hz  Nyquist  frequency. 

Subsequently,  these  five  to  six  power  spectra  were  ensemble  averaged  and  addi¬ 
tionally  smoothed  with  a  20  percent  relative  bandwidth  filter  before  plotting. 


2.4  High-resolution  time  series  plots  of  one-min 
running  mean  values 

Running  mean  caluclations  of  variance  and  important  fluxes  (heat  flux  vnl  shear 
stress)  were  performed.  In  addition  to  the  covariance  matrices  of  U.  V  X ,T  based 
on  the  entire  ‘run’’  (55,000  scans),  we  also  calculated  the  one-min  running  mean 
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of  the  vertical  variance  { ww )  along  the  wind  stress  component  (uu;)  and  the  heat 
flux  {wT).  We  used  a  "box-car’’  filter  function,  i.e.  we  calculated  (ww).  (uw)  and 
( wt )  at  each  time  increment  (0.2  sec',  based  on  the  last  300  scans  (one-min  worth 
of  data)  and  plotted  also  the  result  versus  time  in  a  high-resolution  time  series  plot 
that  can  be  compared  with  the  aligned  measurements  of  V,  V,  IV  and  T.  These 
plots  are  extremely  valuable  for  the  interpretation  of  the  bulk  fluxes  measured. 
They  show  whether  a  measured  10-min  or  92-min  bulk  flux  can  be  attributed  to  a 
single  even*,,  a  few  events,  or  whether  the  measured  quantity  is  evenly  dstributed 
over  the  experiment. 


3  Presentation  of  data 

Next  the  results  of  the  individual  experiments  are  presented.  Each  experiment  is 
represented  by  two  tables  and  five  figures  in  the  following  way. 

3.1  Tables  and  figure  legends 

First  page: 

Log  book  table  for  run  #  N 

The  log  book  table  shows  date  of  experiment,  start  and  stop  time  for  the  KM 
tape  recorder  and  time  period  selected  for  times  enes  and  spectral  analys 

Second  page 

10-min  mean  values  for  run  #  N 

Figures  of  the  variation  of  the  10-min  averaged  surface  layer  scaling  param¬ 
eters  of 

L.  mean  wind  speed 

2.  mean  direction 

3.  mean  temperature 

4.  total  energy 

5.  5 tress 

6.  heat  flux 

7.  stability  parameter  (zj L). 

These  figures  have  already  been  presented  in  part  l  but  here  they  are  repro¬ 
duced  with  the  selected  periods  of  time  series  and  spectral  analysis  marked 
(shaded  periods). 

Third  page 

Sonic  time  series 

Some  55,000  values  of  the  10  Hz  digitized  wind  components  [U,V  W)  are 
shown  after  calibration  and  alignment  ((V')  =  (W)  =  0  0),  together  with  the 
temperature  signal  T.  The  scales  are  U  and  V  :  ±10  m/s;  W  :  ±5  m/s.  and 
TO  —  10 K.  The  mean  temperature  (T)  is  arbitrarily  offsett  from  plot  to  plot 
(in  order  to  obtain  full  resolution  during  recording  and  digitizing). 
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Fourth  page 


Mean  statistics  for  time  senes 

The  table  shows  the  mean  and  covariance  statistics  calculated  for  the  entire 
diffusion  tet  period,  with  a  typical  duration  of  92  minutes  (55.000  scans),  to 
be  considered  as  an  alternative  '‘overall”  average  of  the  10-min  by  10-min 
statistics  for  the  same  period  preented  in  part  I  (see  "second  page’’  above). 

Fifth  page 

Wind  speed  ({/-component)  and  one-min  running  mean  statistics  of  vertical 
variance  (ww),  shear  stress  (uu/).  and  heat  flux  (wT) 

The  first  time  series  plot  out  of  four  is  a  reproduction  of  the  //-component 
(see  “third”  page  above),  but  augmented  here  in  scale  (-2  m/s  to  +8  m/s). 
It  has  been  included  again  for  "event"  marking  and  interpretation. 

The  second  trace  shows  the  one-min  running  mean  value  of  the  vertical  ve¬ 
locity  variance  (ww).  The  scale  is  from  0.0  to  1.0  [mm/ssj. 

The  third  trace  is  the  one-min  running  mean  value  of  the  stress  component 
(uu/).  The  scale  is  from  -1.00  to  +1  00  [mm/ssj. 

The  fourth  and  last  trace  shows  the  one-min  running  mean  value  of  the 
vertical  heat  flux  (uT).  The  scale  is  from  -0.5  to  +0.5  [m/s  K).  When 
multiplying  by  approximately  1200  {p.CP),  this  figure  shows  the  heat  flux 
in  Watts/nr  in  the  interval  between  -(300  and  -H>00  [Watt / nr j. 

Sixth  page 

Velocity  (».  v  and  u)  spectra  for  run  #  N 

The  figure  shows  the  three  velocity  spectra  calculated  as  described  above. 
The  abscissa  ranges  from  0.001  to  10  Hz  in  frequency,  the  Nyquist  frequency 
is  5  Hz.  The  ordinate  shows  the  spectral  density  fS{f)  in  the  interval  be¬ 
tween  0.0001  to  1.0  [mm/ssj.  The  lowest  spectral  estimate  is  at  1/819.2  Hz. 
cf  above. 

Seventh  and  last  page 

Spectrum  of  fluctuating  temperature  for  run  #  N 

The  figure  shows  the  temperature  spectrum  in  the  same  as  for  the  previous 
figure  (velocity  spectra).  Here,  however,  the  ordinate  has  the  dimension  [h':] 
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4  Results 

4.1  Run  #  1,  23  September,  overview 


i  Speed: 
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Figure  I:  10-mm  mean  values  for  Run  #  1. 
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Fi ptrr  i  u.t  unit  H.'- •>pc.:tra  for  Run  pf  l 
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Project  Wire  Pun!  (Very  Convective)  . erroerature  spectrum  (wnae  run) 


Figure.  5:  Temperature,  spectrum  for  Run  if  l. 
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4.2  Run  #  2,  24  September,  Overview 


Speed: 

2.38  cm/s 

(15/16) 

(max  8  hrsj 

ch  l 

ch  2 

ch  3 

ch  5 

AMADEUS  1987  Sonic  Spectra 
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Start  Stop  Duration 


Spectra 


Table  3.  Mean  statistics  for  time  series 


Statistics  from  55200  samples 


Mean 

u  :  0.279  v  : 

0-000 

w  : 

0.000 

C  'ovariance 

,  uu  0.13239  uv  : 

00 10:54 

uw  : 

-0.05721 

uT  : 

0.00050 

vv  : 

0.10501 

vw  : 

-0.00519 

kT 

0.00889 

ww  . 

0.03:157 

wT  : 

-0  02077 

TT  : 

0.08506 
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Figure  S:  Wind  speed  (u)  and  l- 


{ ww ) .  shear  stress  (uw).  and  (sens, hie)  heat  flat  ( wt ) 


min  running  mean  statistics  of  vertical  variance 
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Figurt  0:  u,v  and  w-spcctra  for  Run  tf 
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4.4  Run  #  4,  26  September,  Overview 
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Figure.  11:  10-mm  mean  values  for  Run  #  j. 
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Table  4-  Mean  statistics  for  time  series 


Statistics  from  55000  samples 


j  l 

i  Mean  '  u  :  3.511  v: 

-0.000 

W  ‘ 

-0.000 

Covariance  j  uu  2.72007  uv  : 

0.5401(3 

HW  . 

-0.28901 

uT 

-0.01009  ; 

!  lV  : 

3.3192!) 

IW  : 

0.134(54 

l-T  ; 

0. 33(308  ; 

ww  : 

0.40845 

wT  : 

O.i  3848 

TT 

0.51484  j 
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Figure  15  Temperature  spectrum  for  Run  # 
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Run  #  5,  27  September,  Overview 
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Figure  16:  10- min  values  for  Run  #  5. 
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Table  5.  Mean  statistics  for  time  senes 


Statistics  from  55200  samples 


Mean  j  «  :  0.7-1 1  v 
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Figure  I'J:  u.  r 


>md  w- spectra  for  Run  #  5. 
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figure  JO:  femperature  spec trum  for  Run  # 


i 


'.'Mfil 


4.6 


Run  #  6,  28  September,  Overview 
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Table  ti.  Mean  statistics  for  time  series 
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Run  -#■  7,  30  September,  Overview 
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Figure.  26:  10-min  mean  values  for  Run  #  7. 


16 


Riso-M-2861 


I 


} 

figure  27:  Some  time  series. 


Ris«-M-2861 


i 


h 
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Run  #  8,  30  September,  Overview 
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Figure  31:  10-min  mean  values  for  Run  #  V 
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Table  S.  Mean  statistics  Jo r  time  series. 
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Run  #  9,  1  October,  Overview 
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Figure  17:  Sonic  time  s erie.i. 


Rise-  M  -286 1 
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Run  #  11,  2  October,  Overview 
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Run  #  12,  2  October,  Overview 
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Table  12.  Mean  statistics  for  time  series. 
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Table  IS.  Mean  statistics  for  time  series. 
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Figure  7/  10-min  mean  raluei  for  Run  #  If. 
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Table.  1 4.  Mean  statistics  for  time  senes. 
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Figure  73:  Wind  speed  fu)  and  t-mm  running  mean  statistics  of  vertical  variance 
(wiv).  shear  stress  (ui u),  and  (sensible)  heat  flux  (wt) 
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Figure  ?.5:  Temperature  spectrum  for  Run  #  IJ. 
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